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ABSTRACT 

We present a statistical study of the luminosity functions of galaxies surrounding luminous red 
galaxies (LRGs) at average redshifts (z) = 0.34 and (z) = 0.65. The luminosity functions are derived 
by extracting source photometry around more than 40,000 LRGs and subtracting foreground and 
background contamination using randomly selected control fields. We show that at both studied 
redshifts the average luminosity functions of the LRGs and their satellite galaxies are poorly fitted by 
a Schechter function due to a luminosity gap between the centrals and their most luminous satellites. 
We utilize a two-component fit of a Schechter function plus a log-normal distribution to demonstrate 
that LRGs are typically brighter than their most luminous satellite by roughly 1.3 magnitudes. This 
luminosity gap implies that interactions within LRG environments are typically restricted to minor 
mergers with mass ratios of 1:4 or lower. The luminosity functions further imply that roughly 35% of 
the mass in the environment is locked in the LRG itself, supporting the idea that mass growth through 
major mergers within the environment is unlikely. Lastly, we show that the luminosity gap may be at 
least partially explained by the selection of LRGs as the gap can be reproduced by sparsely sampling 
a Schechter function. In that case LRGs may represent only a small fraction of central galaxies in 
similar mass halos. 

Subject headings: Galaxies: groups: general — Galaxies: evolution — Galaxies: interactions - 
Galaxies: elliptical and lenticular, cD 



1. INTRODUCTION 



Massive galaxies in the nearby universe typically have 
very little cold gas, they host old stellar populations 
and e xhibit extremely low specifi c star formation rates 
(e.g.. iRtberl Il973t IPeletierl Il98l IWorthev et al.l IT992t 



Jrgensen 119991: iTrager et al.l 120001: iKauffmann et^l 



2003HBalogh" et a. Jl2004l:IHogg et al 1 120041: iThomas et al 



20051) . Therefore, studies of these galaxies typically 



find that essentially all of the stellar mass growth takes 
place through mergers and other gravitational interac- 
tions, with the relative importance of each process still 
debatable. For example, while some authors find that 
major dry mergers contribute signifi cantly to the mass 



evolu t ion of mas s ive g a laxies (e.g.. Ivan Dokku m et al 
1999 1 : iTran et al.l 120051: Ivan Dokkuml 120051: iBell et al 
20061 iBovlan-Kolchin et al.1 120061: INaab et al.1 120061: 
Mcintosh et al.l 120081 : iBundv et al.l |2009|), others find 
only a mild con t ribution or none at all (e.g 



Patton et al.l 120021: IBundv et al.l 12006; M asiedi et al 
20061: iWake et all 120061 120081 iMasiedi et al.1 1200 
Tojeiro et al.l 12011). Other studies argue that minor 
mergers and low mass accretion events contribute at 
least some of the stellar mass growth in massive galax- 

Kormendv fe Diorgovskil 
3 



ass g 

ies oy e r a longer timescale ( e . g.Jl 
119891: ISchweizer fe Seitzerl I1992T 



Ivan Dokkuml 2005; 



1 Yale University Astronomy Department, P.O. Box 208101, 
New Haven, CT 06520-8101 USA 

2 Department of Astronomy and Astrophysics, The Pennsyl- 
vania State University, 525 Davey Laboratory, University Park, 
PA 16802, USA 

3 Apache Point Observatory, Apache Point Road, P.O. Box 
59, Sunspot, NM 88349, USA 

4 Center for Cosmology and Particle Physics, New York Uni- 
versity, New York, NY 10003, USA 



Naab et al.l 120071: Bournau d et al.1 12 007: Ka yirai et al 
20081; INaab et all 120091: iBezanson et al.l 120091: iTal et al 



200a iRamos Almeida et al.l 120111: ITal fc van Dokkum 



2011f h The unifying challenge of many of these analy- 



ses is that observations of ongoing mergers typically do 
not recover information regarding the mass ratio, and 
therefore implied growth, of the progenitor galaxies. 

As an alternative to studying galaxy mergers many 
authors have turned to examining the environments in 
which these objects reside. For example, pair counts 
and clustering studies observe the progenitor galaxies 
prior to their merging. By studying the neighbors of 
massive galaxies one in effect probes the reservoir of ob- 
jects with which mergers, and consequently mass growth, 
can occur. However, in studies of groups and clusters 
this technique typically relies on assigning environment 
membership to individual satellites using their inferred 
dista nces from the central galaxy (e.g.. iRamella et al.1 
20001: iChristlein fc Zabludoffl 120031: iMartini et al.1 1200€ ' 



Muzzin et aLM2 009; Wil son et al.ll2009t7chiboucas et al.1 
2010tlTanaka et al.ll2010D . This method requires that sig- 



nificant telescope time is devoted to spectroscopic mea- 
surements of all member candidates in order to properly 
assess their line-of-sight velocity and distance from the 
central galaxy. As a consequence, environments that are 
studied in this way are often restricted to the low redshift 
universe or to small statistical samples for which suffi- 
cient ly de ep data are avai l able (such as the ICoil et all 
120061 and IBolzonella et al.l 120101 studies) . Other tech- 
niques, such as the widely used friends-of-friends algo- 
rithms introduce another complication by explicitly as- 
suming the size, and therefore inferred mass, of candi- 
date groups (most of these algorithms are essentially a 
variation of the technique proposed and demonstrated 



2 



Tal et al. 



m 1e4 
v 

1 le3 

CP 

z 1e2 
1.0 

0.8 

w 

0) 

i; 0.6 
0.4 




0.2 



• 0.28<z<0.40 

• 0.60<z<0.70 




-22 



-23 



g.rest 



-24 1e3 1e4 
N galaxies 



Figure 1. Rest-frame color and absolute magnitude distributions 
of the two redshift samples. The grayscale data points show the 
full SDSS/BOSS LRG data set within the studied redshift range 
while the over-plotted red and blue points represent the density of 
selected galaxies in the high and low redshift bins. The top panel 
shows the initial selection of the most luminous galaxies at high 
redshift and the number density matched samples at low redshift. 
The right panel shows the color matching selection which excluded 
all galaxies that are bluer than the distribution peak value by more 
than lcr. 

in IHuohra fc OeTlerl fl982l and IGeller fc Huchral HQS! . 
Newer algorithms attempt to overcome the need for this 
assumption by iteratively measuring the light distribu- 
tion within a given radius and adjusti ng it using the in- 
ferred mass from a halo density model (|Yang et al.ll2005l 
l2008h . 

Here we follow a different approach and detect satellite 
galaxies in a purely statistical manner, in practice mea- 
suring the average local overdensity of galaxies around 
massive red galaxies. This technique has been applied 
in studies of cluster ga laxies by many au t hors utilizing 
various data sets ( e.g.. IDriver et all 119941 ; IGaidosI 11997 
Lobo et all 119971: ILumsden et all 119971: IValotto et al.| 
1997 1 : I Wilson et al.|[19"97t iFairlev et alj|2002t IGoto et al.l 
200a |Lin et all 120041: iWake et al J 12003: lLoh et al.l 12008. " 
Lu et al.l I2009H . The main advantage of any statistical 



study of this nature is that it requires a minimal set 
of a-priori assumptions regarding the properties of the 
studied objects. Its main disadvantage is that no de- 
tailed information can be extracted regarding any single 
satellite. 

Throughout the paper we adopt the following cosmo- 
logical parameters: Sl m — 0.3, ti\ = 0.7 and H = 70 
km s _1 Mpc -1 . 

2. DATA 

2.1. Sample selection 

We selected galaxy images for this study from the 
seve nth data release of the Sloan Digital Sky Survey 
(see lYork etaLl I2000t IPukugita et all 119961: iGunn et al 



1998; Smith et al 



Gunn et al 



20021: IPier et al.1 12001 _ 
2006; Ab azaiian et al.l l 2009f) and from the interim cata- 
log of the SDSS-III Baryon Oscillation Spectroscopic Sur 



vejO fBOSS. ISdiiegel et al.ll2009t lEisenstein et al.ll201lD . 
This version of the catalog covers an area of roughly 
1620 deg 2 on the sky (compared to 8200 deg 2 in SDSS 
DR7) and it consists of more than 150,000 LRG candi- 
dates with spectra at 0.4 < z < 0.7. All selected ob- 
jects were initially identified by the SDSS pipeline as 
luminous red galaxies from their central surface bright- 
ness and locati on on a rotated color-color d iagram (for 
full details see lEisenstein et al.l l200ll 120 111) . Roughly 
90% of all LRGs are central halo galaxies, thought to 
be residing in gro ups with a typical halo mass of a 
few times 1O 13 M fWake et al.|[200l IZheng et all 120091: 
IReid fc Spergel|[2009|) . The full data set consists of galax- 
ies in two redshift bins as determined spectroscopically 
in the surveys. The low redshift sample consists of SDSS 
LRGs in the redshift range 0.28<z<0.4. These are the 
reddest, most luminous galaxies in the nearby universe 
and they occupy the high mass end of the stellar mass 
spectrum between 1O 11 M and a few times 1O 12 M . The 
high redshift sample is comprised of BOSS "CMASS" 
galaxies, a selection aimed at fin ding objects at an ap- 
proximately const ant stellar mass (|Eisenstein et al.ll2011t 
White et al. 20TTI). The color distribution of these galax- 
ies is broader and this sample includes objects that are 
bluer than the low redshift LRGs. We address this in the 
following subsection. 

2.2. Number density and color matching 

We constructed matched samples at z ~ 0.34 and 
z ~ 0.65 by requiring that galaxies in each redshift bin 
have the same number density and a similar color dis- 
tribution (with an offset due to passive evolution). This 
selection allows us to follow the mass evolution of mas- 
sive galaxies, assuming that the overall number density 
of these objects does not change significantly in the stud- 
ied redshift range. The benefits of a number-density 
limited, rather than luminos ity- or mass-limited , samp le 
have been discussed by e.g., Ivan Dokkum etall (|2O10l) . 

We started by K-correcting the SDSS extinction- 
corrected model magnitudes in the g— and r— bands and 
the BOSS magnitudes in the r— and bands to z = 
following the filter tr a nsform ation technique given by 
Ivan Dokkum fc Franxl (|1996f) . This method finds the 
best-fit linear combination of galaxy templates to an 
observed color and calculates the transformation coef- 
ficients for a given rest-frame filter. We used the peak of 
the observed color distribution as input (g—r = +1.7 and 
r — i = +1.2 for the low and high z samples, respectively) 
and derived rest-frame g— magnitudes and g — r colors for 
all galaxies in the two samples. The imaging bands used 
in this study were selected to minimize the necessary K- 
corrections as the r— band at z = 0.34 and the i— band 
at z = 0.65 roughly overlap with the restframe g— band. 
To test this, we compared the difference in derived K- 
correction values between the two redshifts with stellar 
population synthesis models (|Bruzual &: Chariot! 120031 : 
IMaraston et al.l [20091 ) and found agreement to within a 
few hundredths of a magnitude. 

Our initial selection included objects from the high- 
rcdshift sample that are brighter than M g ^ est — —22.81 
mag. This selection minimizes the effects of incomplete- 
ness due to the BOSS detection limit as it excludes 

5 As of January 20th 2011 
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Figure 2. Examples of LRG and random aperture selections (solid and dashed lines, respectively) for both rcdshift ranges, overlaid on 
a full SDSS imaging field. Photometry is extracted for all objects (including LRGs) in the selected regions within 1000 kpc apertures 
using a threshold of lcr above the background level. Objects are marked with red circles and include LRG satellites as well as foreground 
stars, foreground and background galaxies and spurious detections of noise peaks. We perform 14 and 29 random region selections and 
photometry measurements for the low- and high-redshift LRGs, respectively, to cover a similar area as the full imaging field. 



galaxies close to or fainter than the turnover point in 
N(M g ,rest)- We then excluded the faintest galaxies from 
the low redshift sample until their overall number den- 
sities matched that of the high-z sample to within 1%. 
Lastly, we fit a Gaussian curve to the g — r color distribu- 
tion of the resulting samples, redward of the distribution 
peak, and excluded all the galaxies that are bluer than 
the peak by more than la. Put differently, we selected 
an identical fraction of the reddest galaxies from each 
sample, assuming a normal distribution for both. How- 
ever, since this color selection slightly changed the overall 
number density of the high redshift sample, we iterated 
over the last two steps until both criteria were satisfied, 
resulting in a matched galaxy number density of 4 x 10 -5 
Mpc~ 3 . 

The initial selection criteria of the BOSS "CMASS" 
sources aim at finding galaxies at a roughly constant 
mass. This implies that the BOSS pipeline selects galax- 
ies with a wider range of properties compared to the 
SDSS pip eline, including ga l axies of late type mo rphol- 
ogy (e.g.. IWhite et all l201ll iMasters et aljboilh . The 
color and luminosity that we apply on both BOSS and 
SDSS data sets minimize the contribution of such sources 



to our samples as the final selection only includes the 
reddest, most luminous galaxies at each redshift. We 
therefore refer to these galaxies as "LRGs" throughout 
the paper. We note that some contamination of dusty 
galaxies may still exist in the BOSS data set as simple 
color cuts likely do not complete ly exclude them (e.g., 
IBell et al.ll20nllLabb et al.ll200l . 

The color-magnitude distribution of the initial data 
sets can be seen in Figure [1] along with the selection 
cutoff values and resulting samples. The difference in 
cutoff rcstframe M g between the two redshift samples 
is less than 0.4 magnitudes and roughly corresponds to 
passive luminosity evolution. The final data set includes 
12,813 galaxies at 0.6 < z < 0.7 and 29,477 galaxies at 
0.28<z<0.4. 

3. STATISTICAL DERIVATION OF THE LRG SATELLITE 
LUMINOSITY FUNCTION 

3.1. Region selection and object photometry 

Each imag ing field in SPSS is 8 11" x 590" in angu- 
lar size (see iStoughton et al.l 120021 for a description of 
the SDSS image processing), corresponding to a phys- 
ical scale of roughly 5.6x4.1 Mpc at z = 0.65 and 
3.9x2.8 Mpc at z = 0.34. We acquired r- and 
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i— band images from the survey archives for the se- 
lected galaxies in the low and high redshift samples, re- 
spectivel y, and extracted source photometry using SEx- 
tractor (jBertin fc Arnoutsl [l996h in a 1000 kpc annu- 
lus around the LRGs (145" at z = 0.65 and 207" at 
z = 0.34). We utilized a low detection threshold of la- 
above the background, as determined by SExtractor, re- 
sulting in a large number of spurious detections of nise 
peaks in the images. However, since the random field cat- 
alogs contain the same false detections, this contribution 
to the luminosity function will be taken into account. We 
then calibrated the luminosity of each object using the K- 
correction value and luminosity distance of its adjacent 
LRG, assuming that all sources in any given field are at 
the same redshift. Although this assumption is grossly 
incorrect for any object that is not physically associated 
with the LRG itself, the large numbers of studied sources 
around the LRGs and in random fields ensure that each 
luminosity bin is well sampled in both data sets. We dis- 
cuss this further in Section [3~2l Finally, we divided the 
resulting all-inclusive catalogs into luminosity bins of size 
0.01 dex, and summed the number of objects in each bin. 
A typical region selection at each redshift range can be 
seen in Figure [5J 

3.2. Background and foreground subtraction 

Statistical subtraction of background and foreground 
sources in wide-field surveys can be performed locally 
or globally, with both methods providing successful and 
not inconsistent results (some discussion on the differ- 
ent approaches can be found in lGoto et al.l [20031 and in 
lLoh et all I2008D . We chose to use a local estimate of 
this contribution using source extraction and photometry 
in randomly-selected regions taken from the same SDSS 
fields in which the selected LRGs reside. This method 
better samples the large scale structure in which massive 
LRGs typically reside and therefore attempts to account 
for resulting line-of-sight overdensities. 

For each selected LRG we performed object photome- 
try in a randomly positioned circular region, using iden- 
tical thresholds as those used for the LRG fields (Section 
13. 1|) . The large number of LRGs is well suited to this 
technique as it ensures sufficient sampling of foreground 
and background objects. We allowed the random aper- 
ture centers to lie anywhere in the images such that the 
both LRG and random fields suffer from identical edge 
effects. This ensures that the cumulative size of all ran- 
dom apertures is the same as the total LRG field size. 

We then incorporated the measurements from all fields 
and binned them similarly to the LRG centered regions 
to derive a cumulative source light distribution of fore- 
ground and background sources. We repeated the last 
two steps until we had 10 random iterations for each 
LRG that did not overlap with the LRG aperture itself. 
We also allowed the random apertures to overlap with 
each other by no more than 30% of the apertue area to 
ensure that the errors are not significantly correlated. 

3.3. Photometric errors 

Photometric measurements of the faintest detected ob- 
jects may suffer from significant errors due to statistical 
variation in the background and other residual system- 
atic biases (e.g., flat fielding issues, incorrect sky level 



determination). In order to estim ate these er r ors w e fol- 
lowed the technique suggested bv lLabb et al.l (|2003l ) and 
utilized circular apertures to measure the total flux dis- 
tribution in empty regions in both data sets. First we 
constructed a catalog of all the elliptical apertures that 
had been used by SExtractor to measure photometry 
from the images. We then calculated the average area 
of these apertures in each luminosity bin and produced 
a circular aperture of equivalent area. Next we masked 
out all the objects in each imaging field and measured the 
total flux in every aperture when placed at a random po- 
sition in the field. We repeated the last step 1000 times 
and calculated the standard deviation of the resultant 
measurement distribution for each luminosity bin. 

3.4. Luminosity functions 

The galaxy luminosity functions that were derived in 
Section [3TT] are dominated by light from background and 
foreground sources in the LRG fields. In order to prop- 
erly remove this contamination we averaged the random 
aperture measurements that were carried out in Section 
13. 2\ in each luminosity bin, and subtracted them from 
the LRG measurements. We then normalized the resul- 
tant galaxy luminosity function by the total number of 
groups and clusters, assuming that every studied LRG 
is the central galaxy in its environment. A potential 
concern is that we may preferentially miss LRGs that 
are paired with other LRGs due to incompleteness of the 
SDSS redshift catalogs resulting from fiber collisions. We 
note that such pairs are rare, making up only a few per- 
cent of the sample at the selected lu minosity range, after 
correcting for fiber collisions (see e.g.lMasiedi et al.ll2006t 
Ivan den Bosch et al.ll200l iWake et al.ll2008l ). Therefore, 
the effect on the selection is small. Importantly, all LRGs 
that do not have a redshift measurement are included in 
our analysis as they are in the photometric catalogs that 
we used to construct the luminosity functions. 

In addition, we calculated the difference between SDSS 
model and SExtractor A UTO magnitudes and derived 
a correction factor for the latter. We applied this fac- 
tor to the satellite functions assuming that SDSS model 
magnitudes better assess the total flux from satellite 
galaxies. We further discuss this in the appendix. Fi- 
nally, we applied the same filter transformations dis- 
cussed in Section ^. 2l to the luminosity functions and nor- 
malized the resultant measurements t o solar luminosities 
using Mq g = 5.12 (in the AB system. fBlanton &: Roweisl 
l2lffll . 

The calibrated LRG satellite luminosity functions are 
shown in Figure |3J The blue dot-dashed and the pur- 
ple dashed lines show the observed luminosity functions 
of the LRG and random fields, respectively. The red 
data points show the difference between the two curves 
which is henceforth assumed to represent the luminosity 
function of galaxies associated with LRG environments. 
These luminosity functions also include the central LRGs 
themselves, for which photometry was extracted in an 
identical way to all other sources in the field. The verti- 
cal error bars represent la of the distribution of random 
measurements in each luminosity bin (Section 13. 2|) and 
the horizontal error bars show the standard deviation of 
1000 empty aperture measurements (Section [33)1 . 

The luminosity functions are markedly different from 
a Schechter function and show a clear peak at the lumi- 
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Figure 3. Galaxy luminosity functions of LRG satellites at average redshifts z ~ 0.34 (left panel) and z ~ 0.65 (right panel). The blue and 
purple curves represent the "luminosity" function of the LRG and random fields, respectively, and the red data points show the difference 
between the two. The black curve is a functional fit to the data at luminosities brighter than the detection threshold while the pale orange 
points show data at fainter luminosities. We note that the gap between the LRG and the most massive satellite galaxy is evident in both 
redshift bins. 



nosity of the LRG at both redshifts. This peak is even 
visible in the observed luminosity function (the blue dot- 
dashed curve which is not corrected for foreground and 
background objects). The fact that a large fraction of 
the luminosity in a typical LRG environment is locked 
in the central galaxy itself is the key result of this study. 
In the following sections we will quantify the luminosity 
functions and the total luminosity in the LRG and in its 
satellite galaxies. 

4. QUANTIFYING THE LUMINOSITY FUNCTION OF LRG 
ENVIRONMENTS 

4.1. Functional fits 

The derived luminosity functions that are presented in 
Figure [3] are very different from a Schechter distribution, 
in agreement with halo-b ased studies of nearby and dis- 
tant environments (e.g.. lYang et al.| 120081; IBrown et al.l 
120081) . Following lYang et al.l (|2008l ) we instead quan- 
tified the galaxy luminosity distributions using a two- 
component fit: 



distribution we used a Schechter (1976) function: 



$sat(L) 





L ' 


UJ exp 





(3) 



$(L) = $ LiiG (L)+$ sat (L), 



(1) 



where §lrg and $ sa t are the distributions of LRG and 
satellite luminosities. For the LRG luminosities we as- 
sumed a log-normal distribution, 



LRG 



(L) 



2ircr, 



exp 



(logL - log L c y 
2a? 



(2) 



where the number of LRGs per group or cluster A is set 
to 1 and the distribution width a c and peak center log L c 
are left as free parameters. For the satellite luminosity 



Attempts at fitting the satellite luminosity distribution 
with a modified Schechter function, as suggested by 
lYang et ail (|2008D , resulted in worse agreement with the 
data. We fit all five free parameters (cr c , log L c , <j> a , log L s 
and a s ) simultaneously using t he non-linear least squares 
curve fitting program MPFIT (|MarkwardtJl2009| ). Since 
our data quickly become incomplete below some flux 
threshold we restricted the functional fits to luminosities 
higher than the turnover luminosity of log-L m .; n = 10.1 
and 10.4 in the low and high redshift bins. The best-fit 
parameters are given in Table Q] and the resulting lumi- 
nosity function curves are shown in Figure[3]as black solid 
lines. The slope of the power-law part of the Schechter 
function, represented by the parameter a s , is not well 
constrained and is sensitive to the adopted value of L m i n . 
When varying \og(L m i n ) from 10.1 to 10.6 for the low- 
redshift sample the slope changes by 0.4 dex. To better 
estimate the faint end of the Schechter function we re- 
peated the same model fitting for LRGs in Stripe 82. 

Stripe 82 is a narrow region in the sky that was imaged 
multiple times as par t of SPSS under a varie ty of observ- 
ing conditions (e.g.. [Abazaiian et al.l [2009h . Individual 
images of the same fields were summed to create deep 
co-added frames of this region, each made of 11 to 75 
exposures (with a median of 52 exposures). We repeated 
the analysis described in Section [3] in fields containing 
all 1415 Stripe 82 LRGs at 0.28 < z < 0.40 (redshifts 
were measured as part of the main SDSS survey). For 
the functional fits we kept the parameter a c fixed since 
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Table 1 

Best-fit parameters of the galaxy luminosity function. 





SDSS 


BOSS 


Stripe 82 




0.28 < z < 0.40 


0.60 < z < 0.70 


0.28 < z < 0.40 


A (fixed) 


1.00 


1.00 


1.00 


cr c 


0.09 


0.12 


0.09 (fixed) 


logL c 


11.04 


11.14 


11.00 


<ps 


18.6 


15.3 


24.3 


log L a 


10.5 


10.5 


10.4 


as 


-l.ll±0.47 a 


-0.46±1.15 


-0.95±0.18 



a Errors were calculated by varying the analysis threshold 
log(L min ) by 0.5 dex. 



10.0 



0.1 



1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 


i i i 1 i i 














M \ 






\ • 


0.28 < z < 0.40 




SDSS LRGs (28016 galaxies) 

Stripe 82 LRGs (1415 galaxies) 

. I , , , , I , . , , I , , , , I . , , , I ,' 


* 

1 • s "■ • " 
r •* *." • ■ 

"1 «% 

, 1, , .' 


9.0 9.5 10.0 10.5 11.0 
log(L/L s ) 


11.5 



Figure 4. Comparison between the luminosity functions derived 
from individual SDSS LRG frames (blue data points) and from 
deep Stripe 82 stacks (green data points) in the redshift range 
0.28 < z < 0.40. Solid lines are functional fits to the data using 
the two-component model described in Section \3\ The faint-end 
slope of the Schechter function can be reliably measured and it has 
a value of —0.95. 



the relatively small number of LRGs in Stripe 82 implies 
that the gap region of the luminosity function is not well 
constrained. To calibrate the photometry we matched 
the resulting luminosity of the LRGs themselves to that 
of the same LRGs in the SDSS frames. Figure @] shows 
a comparison between the luminosity function that was 
derived from individual SDSS frames (blue points) and 
the one derived from the deep Stripe 82 stacks (green 
points). As expected, there is good agreement between 
the two curves down to the SDSS threshold luminosity 
of logL ~ 10. ILq. The faint-end slope of the Schechter 
function is better constrained from the deep Stripe 82 
data and it has a value of —0.95. 

4.2. Measurements of the gap width 

The most outstanding feature in both galaxy luminos- 
ity functions that are presented in Figure [3] is a gap at 
the bright end between the LRG luminosity and that of 



the most luminous satellites. Similar luminosity gaps in 
nearby (z < 0.2) massive groups and clusters are typi- 
cally interpreted as a proxy of the magnitude difference 
between t he first and second most luminous grou p mem- 
bers (e.g., IPonman et alj[l99l lYang et al.ll2008[) . How- 
ever, the statistical nature of this study and the methods 
used to derive the galaxy luminosity functions imply that 
we cannot measure such a magnitude difference for any 
given group or cluster as membership is not assigned to 
individual sources. Instead, we treat the gap in the lumi- 
nosity function as a probability distribution for finding 
satellites at a given relative luminosity compared to the 
central. Thus, we quantify the magnitude gap by find- 
ing the luminosity above which LRG groups and clusters 
have on average exactly one satellite. The ratio between 
this luminosity and the peak of LRG luminosities is then 
roughly equivalent to the magnitude gap between the 
two most luminous members of the environment. Put 
differently, we calculate the luminosity L u above which 
the integrand over the satellite luminosity distribution 
equals unity: 



(«., + !) 



— exp 



L 



dlogL^l (4) 



This statistic implies a gap width of log L c — log L u ^0.5 
dex, or roughly 1.3 magnitudes, at both redshifts. An al- 
ternative measure of the gap width is the difference be- 
tween log L c and log L s which is also consistent between 
the two redshift samples and has a value of roughly 0.5 
dex. We note, however, that this measurement is less ro- 
bust as the parameter log L s is degenerate with the other 
parameters of the Schechter function. 

5. DISCUSSION 

5.1. Evidence for an early formation of LRG 
environments 

Measurements of the magnitude difference between 
the first and second most luminous group members have 
typica lly been interprete d as a gauge of the group age 
'e.g., iSandaee fc Hardvl 119731: iTremaine fc Richstone 



1977 : ISchneider et al.lll983t !Barnedll98ft IPonman et al 



2005 



2007 



1994 iKhosroshahi et all 120041: iD'Onghia et all 
Milo savlievi et all 20061: van den Bosch et al.l 
Yang et al . 2008r iDariush" et al J 1201(1 ). In this model 
the most massive members merge quickly (within a 
few tenths of a Hubble time) an d leave behind only 
significantly less massive satellites. iMasiedi et al.l (|2008l ) 
showed supporting evidence for this model by deriving 
the small scale correlation function of SDSS LRGs and 
estimating that these central galaxies can only grow by 
up to a few percent via merging with their satellites. 
The galaxy luminosity functions presented in figure [3] 
further support this and suggest that LRGs typically 
live in groups where the central galaxy is significantly 
more massive than its most luminous satellite. 

In addition, Figure [5] shows that the depth of the lumi- 
nosity gap at z ~ 0.34, as well as its width, are generally 
similar to the gap properties at z ~ 0.65. This is consis- 
tent with no significant evolution in the gap properties 
between the two redshift bins and therefore no signifi- 
cant merger activity between six and four billion years 
ago. Moreover, the existence of the gap at z = 0.67 is 
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Figure 5. Luminosity distributions of the two rcdshift bins over 
plotted with their respective best fit models. The gap in the lumi- 
nosity function is evident, and generally similar in properties, at 
both redshifts. The ratio of LRG pea k luminosities i s cons istent 
with the luminosity evolution model of Blanton ct al. (2003). 

consistent with the scenario in which the gap (and group) 
formed early. 

We note that the ratio between the two LRG peaks 
is consistent with passive lumino sity evolution as de- 
scribed by the model derived by iBlanton et al.l (2003) 
for the same galax i es in SDSS and is in agreement with 
To ieiro fc Percivall (|2010h . 

5.2. Mergers and the mild mass evolution of LRGs 

The luminosity functions derived in Sections |3] and 
[4] can be used to directly characterize the reservoir of 
satellite galaxies with which the central LRG can merge. 
We utilize these luminosity distributions to calculate the 
mass of satellite galaxies which contribute the most to 
the LRG stellar mass growth. In Figure[S]we present the 
distribution of total galaxy luminosities per log luminos- 
ity bin for the low-redshift sample. From this relation 
we can measure the ratio between the peak of LRG and 
satellite galaxy luminosity distributions. This ratio has 
a value 4:1 and it suggests that the satellite mass distri- 
bution peaks at 25% of the LRG mass. This is consistent 
with the derivations of gap width fSection l4.2l) which im- 
ply that the most luminous LRG satellite is on average 
more than three times fainter than the LRG itself. This 
result supports a scenario in which major mergers within 
LRG environments are improbable and that any mass 
growth takes place t hrough mergers with mass ratios 
of 1:4 or lower (e.g., iMasiedi et al.l 120081: iKavirai et al.l 
120091) . 

In addition, by separately integrating the total LRG 
and satellite luminosity distributions of Figure [5] we cal- 
culate that roughly 35% of the mass in LRG environ- 
ments is locked in the central galaxy itself. Although 
the potential for stellar mass growth via mergers within 
groups or clusters is not insignificant, significant evolu- 
tion, when restricted to growth through minor mergers, 
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Figure 6. Total luminosity per log luminosity bin for the low- 
redshift sample. This figure shows that the ratio between the peaks 
of LRG and satellite galaxy luminosity distributions is roughly 4:1. 
This ratio implies that most of the satellite mass is in galaxies that 
are roughly four times less massive than the LRG itself. 

may in fact require many Gyrs. For comparison, in mas- 
sive clusters, where the galaxy luminosity distribution 
is better described by a Schechter function, the central 
galaxy contains as little as 15% or less of the total stel- 
lar mass in the environment. In such environments the 
central galaxy can grow quickly in mass through major 
mergers with its most massive satellites. 

We note that this estimate is an upper limit for the 
maximum mass growth through mergers within the LRG 
group. Since the colors of satellite galaxies are expected 
to be bluer than the those of central LRGs, a more ac- 
curate mass calibration will likely shift the satellite mass 
function to lower masses and increase the luminosity gap 
width. Proper estimates of the color of satellites in LRG 
group are therefore important for this analysis but they 
are unfortunately beyond the scope of this study. 

5.3. Uncertainties in mass estimates 

We note that the luminosity, and therefore mass, cal- 
ibrations of the luminosity functions that are presented 
in Figures [3] and [S] are only as valid as some of the as- 
sumptions that were made for the functional fits. For 
example, in order to calibrate the luminosity functions 
we assumed that LRG satellites typically have the same 
colors as the central galaxy. If the satellites are instead 
bluer than the LRG, their calibrated luminosities would 
be fainter, implying lower total stellar mass. In that case 
the mass function of LRG environments would exhibit an 
even wider gap between the two most massive members. 
Although group satellites have indeed been observed to 
be bluer than a typical LRG, most of them are close 
in color to their central galaxy (e.g.. lYang et al.l 120081 : 
iBalogh et al.ll2009t IWeinmann et alll2009D . 

An additional error may arise from the magnitude cor- 
rection that was derived in Section 13.41 and is presented 
in Figure [4] This estimate relies on the assumption that 
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Figure 7. Simulated galaxy luminosity functions for group and 
cluster sizes of 30,200 and 1000 galaxies per environment. Galaxy 
luminosities were drawn at random from an assumed underlying 
distribution which is modeled as a Schcchtcr function. Each curve 
represents the average of 1000 groups and clusters where at least 
one galaxy is more luminous than a threshold of logL = 11. The 
shape of the observed galaxy luminosity function is highly corre- 
lated with the number of satellites even for a constant underlying 
distribution, resulting from the steepness of this distribution at 
high luminosities. Also plotted is the derived luminosity function 
for the SDSS data set of LRGs at 0.28 < z < 0.40. 

the magnitude difference between the two photometric 
measurements remains constant even at fainter bright- 
nesses. An increase in this correction factor at fainter 
magnitudes may imply that we underestimate the lumi- 
nosity of satellite galaxies and therefore potentially over- 
estimate the width of the luminosity gap. 

Finally, it is possible that sources that are at a small 
angular distance from the LRGs are missed from our cat- 
alogs due to imperfect source separation by SExtractor. 
Although this could also potentially lead to an overesti- 
mate of the gap width it is more likely that sources that 
are missed are faint and do not contribute much to the 
luminosity functions in the gap region. 

5.4. The Schechter function and luminosity gap 

Gaps in the luminosity functions of groups and clus- 
ters are often observed i n studies of the environments 
of lum inous galaxies (e.g.. I Yang et al. 2008; B rown et al.l 
12008ft . These gaps may be created during the initial for- 
mation of massive centrals by mergers of the most mas- 
sive members in an environment. Such (major) mergers 
can increase the luminosity of the central galaxy while 
at the same time effectively decrease the luminosity of 
the most massive satellite. For example, a sparse group 
hosting two L* galaxies would not be selected for this 
study and it would also not have a gap in its luminos- 
ity function. However, some time after the two galaxies 
merge the remnant galaxy will likely be luminous enough 
to be classified as an LRG and a gap in the luminosity 
function will appear. In this scenario, luminous massive 
galaxies in sparse groups are not likely to have any re- 
maining massive companions in their environments. It 



therefore follows that group LRGs preferentially live in 
environments where a luminosity gap can be observed 
even if the underlying luminosity distribution does not 
exhibit a gap. 

Here we test whether the observed gap can be a re- 
sult of the selection of the LRGs, combined with sparse 
sampling of an underlying Schechter function. We start 
by assuming that this distribution is the same in all en- 
vironments, regardless of other properties, and that it 
can be well described by a Schechter function. We then 
randomly sample the assumed Schechter distribution N 
times, where each reading represents a galaxy observa- 
tion and N is the total number of galaxies in a "group" . 
We record the resultant galaxy luminosities and derive 
luminosity functions for all groups. Finally, we repeat 
the last two steps until at least 1000 groups (or clusters) 
of each size N have one or more luminous galaxies with 
log(L) > 11L© and exclude all other groups. Figure [7] 
shows the average luminosity function of the 1000 ran- 
dom samples for group sizes of N = 30, 200 and 1000 
galaxies. The underlying distribution from which galax- 
ies were sampled is also plotted (black dashed line). The 
luminosity gap can be seen in the average distribution of 
groups with 30 and 200 galaxies, but not in the richest 
clusters. 

It is therefore evident that the gap in the luminosity 
function may indeed result from the criteria used to se- 
lect our sample. The steep bright end of the Schechter 
function implies that the probability of having multiple 
luminous galaxies in one environment decreases with the 
number of physically associated members. It then fol- 
lows that LRG environments are not necessarily unique 
in their underlying luminosity distribution but rather in 
that they host a luminous galaxy. Furthermore, this may 
suggest that groups that do not host a bright galaxy can 
easily be missed by group finding studies as all of the 
member galaxies may be undetected. We note that this 
result is not restricted to a Schechter distribution and 
that any underlying distribution which is much steeper 
at its bright end may reproduce a similar luminosity gap. 

Consequently, if environments are selected by their to- 
tal halo mass rather than the luminosity of their most 
massive member, they may exhibit a much shallower gap 
or even none at all. When assuming a universal under- 
lying Schechter distribution for all halos of a given mass, 
groups hosting a massive galaxy become rare and their 
contribution to the average luminosity function dimin- 
ishes. In fact, in our simulations such groups account for 
only 3% of the total environments with the same num- 
ber of satellites. We note, however, that these simula- 
tions are over-simplistic and that they do not represent 
the true global distribution of halo masses. Moreover, 
all groups hosting the same number of galaxies do not 
have the same total mass, suggesting that this value may 
be a gross underestimate of the fraction of groups host- 
ing an LRG. Nevertheless, any selection which is based 
on the central galaxy luminosity would likely underesti- 
mate the scatter in the masses of centrals in s i milar halos. 
Such a scatter was observed by lYang et al.l ((2008), who 
noticed that the distribution of central galaxy luminosi- 
ties (a c ) has a constant width of roughly 0.15 dex for 
Mh > 10 A -1 M©. Moreover, some scatter is required 
to provide good halo model fits to the clustering and 
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space density of LRGs ()Wake et all 120081: iZheng et ail 
l200l IWhite etaLll20ll . 

6. SUMMARY AND CONCLUSIONS 

In this paper we utilized imaging data of more than 
40,000 SDSS LRGs in the redshift ranges 0.28 < z < 0.40 
and 0.60 < z < 0.70 to study the luminosity function of 
LRG satellite galaxies. We extracted source photome- 
try in 1000 kpc apertures around the LRGs themselves, 
as well as in randomly selected fields, to characterize 
and remove background and foreground contamination 
in a statistical way. The large size of SDSS makes it 
an ideal data set for this study as it reduces statistical 
noise and allows for excellent determination of the varia- 
tion between fields. The addition of BOSS spectra to the 
imaging data allowed us to study the galaxy luminosity 
function around LRGs out to z = 0.7 for the first time. 

We successfully determined the galaxy luminosity func- 
tion of LRG satellites down to two magnitudes fainter 
than the LRG brightness at z ~ 0.65 despite the rela- 
tively shallow depth of SDSS. We did so without acquir- 
ing additional spectra to confirm member candidates and 
by assuming that the properties of LRG satellites are 
correlated with the properties of the LRGs themselves. 
This technique can be utilized to study the environment 
around any type of galaxy in a well-defined sample and 
in a sufficiently large data set. 

We found that the luminosity function of LRG envi- 
ronments is markedly different from a Schechter function 
and that LRGs are typically ~ 1.3 magnitudes brighter 
than the next brightest galaxy in their environment. We 
showed that most of the group mass that is not in the 
central itself is concentrated in satellite galaxies that are 
roughly four times less luminous than the LRG. This im- 
plies that major mergers within LRG groups and clusters 
are improbable and that any stellar mass growth likely 
occurs through minor mergers. We also demonstrated 
that the luminosity function gap is already in place at 
z = 0.65 and estimated that the total mass in LRGs ac- 
counts for roughly 35% of the total stellar mass in their 
environment. The existence of the gap at this redshift 
supports that LRG environments typically formed early 
and were already in place 6 Gyrs in the past. 

Lastly, we performed simple assembly simulations of 
galaxy groups and clusters and demonstrated that the se- 
lection criteria used in this study may have preferentially 
picked environments that have a gap in their luminosity 
function. In this case the luminosity gap results from 
the properties of the underlying Schechter function that 
was used in the simulations and from the requirement to 
have a luminous galaxy in every selected environment. 
This further implies that the LRG luminosity function 
may not be inherently unique and that a true halo mass 
selected sample may exhibit a large scatter in the prop- 
erties of its most massive galaxies. 
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APPENDIX 

COMPARISON BETWEEN SDSS AND S EXTRACTOR MAGNITUDES 

Aperture based photometry can lea d to a systematic misestimate of total flux measurements compared to model 
based photometry. I Taylor et al. (2011) performed an extensive comparison of sever al flux estimates an d showed that 
SExtractor AUTO magnitudes are up to 10% fainter than SDSS model magnitudes. I Taylor et al.l (|2011l) further found 
that this value is dependent on the light profile slope of the best-fit galaxy model. 

In order to correct for this effect we extracted photometry for all objects in 150 SDSS fields in Stripe 82 and compared 
the resulting AUTO magnitudes to model magnitudes from the Stripe 82 catalogs. This comparison is shown in 
figure HI where the difference between the measurements is plotted against the SExtractor values. We then fitted a 
line to the running median of this distribution and derived a correction relation for the extracted magnitudes used 
throughout this paper, assuming that Stripe 82 model magnitudes better reflect the total flux of the studied galaxies. 
The correction factor is essentially constant and it has a value of roughly 0.04 magnitudes between msE = 17 and 
m SE — 22. The distribution plo tted in Figure El is c onsistent with the relation between SDSS model and SExtractor 
AUTO measurements found by iTavlor et al.1 (|2011l ) for objects brighter than roughly 17.5 magnitudes. Mill et al~l 
(|2011[ ) compared SDSS and SExtractor flux values over a broad range of magnitudes and also found a small offset 
between the measurements. 




SEKtrocior r-mogniture 



Figure 8. Difference between SExtractor AUTO magnitudes and SDSS model magnitudes for sources in 150 SDSS Stripe 82fields. The 
black data points represent the measured values of all sources while the red thick points are the running median of the sample. The dashed 
red line is a linear fit to the running median and it shows that SExtractor underestimates the total object flux compared to Stripe 82 by 
roughly 0.04 magnitudes. Assuming that SDSS model magnitudes are a better assessment of the total object flux we use this fit to correct 
our SExtractor measurements. 

We note that since all galaxies in our sample do not have the same light profile slope, we introduce a error of up to 
10% to the measured luminosities by applying a single correction relation for all galaxies. This error could potentially 
translate to a systematic offset in the flux measurement of the LRGs themselves if the LRG light profile is significantly 
steeper than the typical satellite profile. However, since the absolute fraction of missed light at any given magnitude 
is small, we expect this contribution to be insignificant compared to other sources of scatter in this study. 



